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The idea that events obey a definite causal order is deeply rooted in our understanding of the world and at 
the basis of the very notion of time. But where does causal order come from, and is it a necessary property 
of nature? We address these questions from the standpoint of quantum mechanics in a new framework for 
multipartite correlations which does not assume a pre-defined global causal structure but only the validity of 
quantum mechanics locally. All known situations that respect causal order, including space-like and time-like 
separated experiments, are captured by this framework in a unified way. Surprisingly, we find correlations that 
cannot be understood in terms of definite causal order. These correlations violate a 'causal inequality' that is 
satisfied by all space-like and time-like correlations. We further show that in a classical limit causal order always 
arises, which suggests that space-time may emerge from a more fundamental structure in a quantum-to-classical 
transition. 



Introduction 

One of the striking features of quantum mechanics is that it 
challenges the view that physical properties are well defined 
prior to and independent of their measurement. This motivates 
an operational approach to the theory, where primitive labora- 
tory procedures, such as measurements and preparations, are 
basic ingredients. Although significant progress has recently 
been made in this direction 1 1 - 8 |, most approaches still retain 
a notion of space-time as a pre-existing 'stage' in which events 
take place. Even the most abstract constructions, in which no 
explicit reference to space-time is made, do assume a definite 
order of events: if a signal is sent from an event A to an event 
B in the run of an experiment, no signal can be sent in the 
opposite direction in that same run. But are space, time, and 
causal order truly fundamental ingredients of nature? Is it pos- 
sible that, in some circumstances, even causal relations would 
be 'uncertain', similarly to the way other physical properties 
of quantum systems are |9 |? 

Here we show that quantum mechanics allows for such a 
possibility. We develop a framework that describes all corre- 
lations that can be observed by two experimenters under the 
assumption that in their local laboratories physics is described 
by the standard quantum formalism, but without assuming that 
the laboratories are embedded in any definite causal struc- 
ture. These include non- signalling correlations arising from 
measurements on a bipartite state, as well as signalling ones, 
which can arise when a system is sent from one laboratory 
to another through a quantum channel. We find that, sur- 
prisingly, more general correlations are possible, which are 
not included in the standard quantum formalism. These cor- 
relations are incompatible with any underlying causal struc- 
ture: they allow performing a task — the violation of a 'causal 
inequality' — which is impossible if events take place in a 
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causal sequence. This is directly analogous to the famous 
violation of local realism: quantum systems allow perform- 
ing a task — the violation of Bell's inequality 1 10] — which is 
impossible if the measured quantities have pre-defined local 
values. The inequality considered here, unlike Bell's, con- 
cerns signalling correlations: it is based on a task that in- 
volves communication between two parties. Nevertheless, 
it cannot be violated if this communication takes place in a 
causal space-time. Previous works about relativistic causal- 
ity in quantum mechanics focused on non- signalling correla- 
tions between space-like separated experiments or on a finite 
speed of signalling |11-19|. In the present work we go be- 
yond such approaches since we do not assume the existence 
of a space-time (or more generally of a definite causal struc- 
ture) on which the evolution of quantum systems and the con- 
straints given by relativity are defined. One of the motivations 
for our approach comes from the problem of time in attempts 
to merge quantum theory and general relativity into a more 
fundamental theory ll20H25ll . 



Results 

Causal inequality 

The general setting that we consider involves a number of 
experimenters — Alice, Bob and others — who reside in sepa- 
rate laboratories. At a given run of the experiment, each of 
them receives a physical system (for instance, a spin-^ par- 
ticle) and performs operations on it (e.g. measurements or 
rotations of the spin), after which she/he sends the system out 
of the laboratory. We assume that during the operations of 
each experimenter, the respective laboratory is isolated from 
the rest of the world — it is only opened for the system to come 
in and to go out, but between these two events it is kept closed. 
It is easy to see that, under this assumption, causal order puts 
a restriction on the way in which the parties can communicate 
during a given run. For instance, imagine that Alice can send a 
signal to Bob. [Formally, sending a signal (or signalling) is the 
existence of statistical correlations between a random variable 
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that can be chosen by the sender and another one observed by 
the receiver.] Since Bob can only receive a signal through the 
system entering his laboratory, this means that Alice must act 
on her system before that. But this implies that Bob cannot 
send a signal to Alice since each party receives a system only 
once. Therefore, bidirectional signalling is forbidden. 

Consider, in particular, the following communication task 
to be performed by two parties, Alice and Bob. After a given 
party receives the system in her/his laboratory, she/he will 
have to toss a coin (or use any other means) to obtain a random 
bit. Denote the bits generated by Alice and Bob in this way by 
a and b, respectively. In addition. Bob will have to generate 
another random bit b\ whose value, or 1, will specify their 
goal: if b' = 0, Bob will have to communicate the bit b to 
Alice, while if b^ = 1, he will have to guess the bit a. With- 
out loss of generality, we will assume that the parties always 
produce a guess, denoted by x and y for Alice and Bob re- 
spectively, for the bit of the other (although the guess may not 
count depending on the value of Their goal is to maximize 
the probability of success 

Psucc •= \ \P{x = bW = 0) + P{y = aW = D] . (1) 

If all events obey causal order, no strategy can allow Alice 
and Bob to exceed the bound 



< 3/4. 



(2) 



Indeed, as argued above, in any particular order of events, 
there can be at most unidirectional signalling between the par- 
ties, which means that at least one of the following must be 
true: Alice cannot signal to Bob, or Bob cannot signal to Al- 
ice. Consider, for example, a case where Bob cannot signal 
to Alice. Then, if b' = 1 , they could in principle achieve up 
to P(y = a\b' = 1) = 1 (for instance, if Alice operates on her 
system before Bob, she could encode information about the 
bit a in the system and send it to him). However, if b^ = 0, the 
best guess that Alice can make is a random one, resulting in 
P(x = b\b' = 0) = 1/2 (see Fig.[T^). Hence, the overall proba- 
bility of success in this case will satisfy psucc ^ 3/4. The same 
holds if Alice cannot signal to Bob. It is easy to see that no 
probabilistic strategy can increase the probability of success. 

Formally, the assumptions behind the causal inequality ^ 
can be summarized as follows: 

Causal structure (CS) — The main events in the task (a sys- 
tem entering Alice's/Bob's laboratory, the parties obtaining 
the bits a, b, and b\ and producing the guesses x and y) are 
localized in a causal structure. [A causal structure (such as 
space-time) is a set of event locations equipped with a partial 
order < that defines the possible directions of signalling. If 
A < B, we say that A is in the causal past of B (or B is in 
the causal future of A). In this case, signalling from A to B 
is possible, but not from B to A. For more details on causal 
structures, see Appendix.] 

Free choice (FC) — Each of the bits a, b, and b' can only be 
correlated with events in its causal future (this concerns only 
events relevant to the task). We assume also that each of them 
takes values or 1 with probability 1/2. 
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Fig. 1 : Strategy for accomplishing communication task by using 
processes with definite and indefinite causal order, (a) There ex- 
ists a global background time according to which Alice's actions are 
strictly before Bob's. She sends her input a to Bob, who can read 
it out at some later time and give his estimate y = a. However, Bob 
cannot send his bit b to Alice as the system passes through her labora- 
tory at some earlier time. Consequently, she can only make a random 
guess of Bob's bit. This results in a probability of success of 3/4. (b) 
If the assumption of a definite order is dropped, it is possible to de- 
vise a resource (i.e. a process matrix W) and a strategy that enables 
a probability of success ^^r^ > 3/4 (see text). 



Closed laboratories (CL) — Alice's guess x can be corre- 
lated with Bob's bit b only if the latter is generated in the 
causal past of the system entering Alice's laboratory. Analo- 
gously, y can be correlated with a only if a is generated in the 
causal past of the system entering Bob's laboratory. 

In the Appendix, we present a formal derivation of the in- 
equality from these assumptions. 

Interestingly, we will see that if the local laboratories are 
described by quantum mechanics, but no assumption about 
a global causal structure is made (Fig. [TJ)), it is in principle 
possible to violate the causal inequality in physical situations 
in which one would have all the reasons to believe that the bits 
are chosen freely and the laboratories are closed. This would 
imply that CS does not hold. 



Framework for local quantum mechanics 

The most studied, almost epitomical, quantum correlations 
are the non- signalling ones, such as those obtained when Al- 
ice and Bob perform measurements on two entangled sys- 
tems. Signalling quantum correlations exist as well, such as 
those arising when Alice operates on a system which is subse- 
quently sent through a quantum channel to Bob who operates 
on it after that. The usual quantum formalism does not con- 
sider more general possibilities, since it does assume a global 
causal structure. Here we want to drop the latter assumption 
while retaining the validity of quantum mechanics locally. For 
this purpose, we consider a multipartite setting of the type out- 
lined earlier, where each party performs an operation on a sys- 
tem passing once through her/his laboratory, but we make no 
assumption about the spatio-temporal location of these exper- 
iments, not even that there exists a space-time or any causal 
structure in which they could be positioned (see Fig. [2]). Our 
framework is thus based on the following central premise: 
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Fig. 2: Local quantum experiments with no assumption of a pre- 
existing background time or global causal structure. While the 
global causal order of events in the two laboratories is not fixed 
in advance and in general not even definite (here illustrated by the 
'shifted' relative orientation of the two laboratories), the two agents, 
Alice and Bob, are each certain about the causal order of events in 
their respective laboratories. 

Local quantum mechanics — The local operations 
of each party are described by quantum mechan- 
ics. 

More specifically, we assume that one party, say Alice, can 
perform all the operations she could perform in a closed labo- 
ratory, as described in the standard space-time formulation of 
quantum mechanics. These are defined as the set of quantum 
instruments |26| with an input Hilbert space *H^^ (the system 
coming in) and an output Hilbert space (the system going 
out). (The set of allowed quantum operations can be used as 
a definition of 'closed quantum laboratory' with no reference 
to a global causal structure.) A quantum instrument can most 
generally be realized by applying a joint unitary transforma- 
tion on the input system plus an ancilla, followed by a projec- 
tive measurement on part of the resulting joint system, which 
leaves the other part as an output. (From the point of view 
of each party, the input/output systems most generally corre- 
spond to two subsystems of the Hilbert space associated with 
the local laboratory, each considered at a diflferent instant — the 
time of entrance and the time of exit, respectively — where the 
subsystems and the respective instants are independent of the 
choice of operation that connects them.) When Alice uses a 
given instrument, she registers one out of a set of possible out- 
comes, labeled by j = 1, . . . , ^. Each outcome induces a spe- 
cific transformation from the input to the output, which corre- 
sponds to a completely positive (CP) trace-nonincreasing map 
EH : £(!H^') £(!H^'\ where XCK^), X = Ai, A2, is 
the space of matrices over a Hilbert space of dimension 
dx. The action of each Mj on any matrix cr e JL^IH^^ ) can be 

written as 1271 M^{cr) = Y!k=i ^jk^^]v ^ ~ <^Ai<^A2. where 
the matrices Ejk : 9^^' "H^^ satisfy S^^i E^^Ejk < t^\Sj. 
If the operation is performed on a quantum state described 
by a density matrix p, M^{p) describes the updated state af- 
ter the outcome j up to normalization, while the probability 
to observe this outcome is given by P{M]) = Tr[Mj(p)]. 
The set of CP maps {At^j _^ corresponding to all the possi- 



ble outcomes of a quantum instrument has the property that 
Mj is CP and trace-preserving (CPTP), or equivalently 

^]k^jk = ' which reflects the fact that the prob- 
ability to observe any of the possible outcomes is unity. A 
CPTP map itself corresponds to an instrument with a single 
outcome which occurs with certainty. 

In the case of more than one party, the set of local outcomes 
corresponds to a set of CP maps Mf, Alf , • • • . A complete 

list of probabilities P (Alf , M^, • • • ) for all possible local out- 
comes will be called process. (It is implicitly assumed that the 
joint probabilities are noncontextual, namely that they are in- 
dependent of any variable concerning the concrete implemen- 
tation of the local CP maps. For example, the probability for 
a pair of maps Mf, to be realized should not depend on 

the particular set {Atf , Atf , Al^} of possible CP maps 
associated with Alice's operation.) A process can be seen as 
an extension of the notion of state as a list of probabilities for 
detection results |3 1 described by a positive operator-valued 
measure (POVM), which takes into account the transforma- 
tion of the system after the measurement and can thus capture 
more general scenarios than just detection. Here we will con- 
sider explicitly only the case of two parties (the generaliza- 
tion to arbitrarily many parties is straightforward). We want 
to characterize the most general probability distributions for a 
pair of outcomes /, j, corresponding to CP maps A(f , A(J, to 
be observed, that is, to characterize all bipartite processes. 

In quantum mechanics, operations obey a specific algebraic 
structure that reflects the operational relations between labo- 
ratory procedures |3|. For example, a probabilistic mixture 
of operations is expressed as a linear convex combination of 
CP maps. It can be shown (see Appendix) that the only prob- 
abilities consistent with the algebraic structure 
of local quantum operations are bilinear functions of the CP 
maps Mf and My Thus the study of the most general bi- 
partite quantum correlations reduces to the study of bilinear 
functions of CP maps. 

It is convenient to represent CP maps by positive semi- 
definite matrices via the Choi-Jamiolkowsky (CJ) isomor- 
phism l^aS. The CJ matrix Mf e £(9^^^ 9^^^) corre- 
sponding to a linear map Mi : JL{9i^^) £(9{^^) is defined 
as Mf^^' := [j0M (ir><n)]\ where = Ijj) ^ 
TY^i (g) is a (not normalized) maximally entangled state, 
the set of states is an orthonormal basis of 9(^\ I 

is the identity map, and T denotes matrix transposition (the 
transposition, absent in the original definition, is introduced 
for later convenience). Using this correspondence, the prob- 
ability for two measurement outcomes can be expressed as a 
bilinear function of the corresponding CJ operators as follows: 

P {Mf, M^) = Tr (Mf Mf ^^)] , (3) 

where \v^i^2fiifi2 ^ matrix in £(9^^' 9^^^ 9(^' 9(^^). 

The matrix W should be such that probabilities are non- 
negative for any pair of CP maps Mf, M^. We require that 
this be true also for measurements in which the system inter- 
acts with any system in the local laboratory, including sys- 
tems entangled with the other laboratory. This implies that 
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W^^^^^^^^ must be positive semidefinite (see Appendix). Fur- 
thermore, the probability for any pair of CPTP maps M^, 
to be realized must be unity (they correspond to instruments 
with a single outcome). Since a map is CPTP if and only 
if its CJ operator satisfies M^'^^ > and Tva^^^'^^ = 
(similarly for M^), we conclude that all bipartite probabilities 
compatible with local quantum mechanics are generated by 
matrices W that satisfy 

^AiA2fiifi2>Q [non-negative probabilities], (4) 

[probabilities sum up to 1]. 

We will refer to a matrix w^^^^^^^^ that satisfies these condi- 
tions as a process matrix. Conditions equivalent to Eqs. ^ 
and ([5]) were first derived as part of the definition of a 'quan- 
tum comb' 1 30 1, an object that formalizes quantum networks. 
Combs, however, are subject to additional conditions fixing a 
definite causal order, which are not assumed here. 

A process matrix can be understood as a generalization of 
a density matrix and Eq. ^ can be seen as a generalization 
of Born's rule. In fact, when the output systems A2, B2 are 
taken to be one-dimensional (i.e. each party performs a mea- 
surement after which the system is discarded), the expression 
above reduces to p(Atf,M^) = Tr [w"^' ^ ' [Mf' <^ M^)], 
where now Mf, are elements of local POVMs and W^^^^ 
is a quantum state. This implies that a quantum state p^i^i 
shared by Alice and Bob is generally represented by the pro- 
cess matrix \y^i^2fiifi2 = pAifii ^ ^MB2 signalling correlations 
can also be expressed in terms of process matrices. For in- 
stance, the situation where Bob is given a state p^^ and his 
output is sent to Alice through a quantum channel C, which 
gives P(Mf, Mp = Tr [Mf o C o (p^i)], is described by 
W"^'"^^^'^^ = ]1^2(g)(c^2Ai)r^pfii ^^j^^j,^^fi2Ai is the CJ matrix 
of the channel C from ^2 to Ai . 

The most general bipartite situation typically encountered 
in quantum mechanics (i.e. one that can be expressed in terms 
of a quantum circuit) is a quantum channel with memory, 
where, say. Bob operates on one part of an entangled state and 
his output plus the other part is transferred to Alice through 
a channel. This is described by a process matrix of the form 
1^2 jyAifiifi2 Conversely, all process matrices of this form 
represent channels with memory |30|. This is the most general 
situation in which signalling from Alice to Bob is not possi- 
ble, a relation that we will denote by A ^ 5 in accord with 
the causal notation introduced earlier. Process matrices of this 
kind will be denoted by W^^^ (note that for non-signalling 
processes, both A^^B and 5 ^ A are true). As argued earlier, 
if all events are localized in a causal structure and Alice and 
Bob perform their experiments inside closed laboratories, at 
most unidirectional signalling between the laboratories is al- 
lowed. In a definite causal structure, it may still be the case 
that the location of each event, and thus the causal relation be- 
tween events, is not known with certainty. A situation where 
B ^ A with probability < ^ < 1 and A ^ B with probability 
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Fig. 3: Terms appearing in a process matrix. A matrix satisfying 
condition (|4| can be expanded as ^^1^2^152 - Yj^vAj ^juv/iycr^^ ®o\^ 

cr^^ (8) cTy'^ , w^yxy ^ where the set of matrices [cr^ ^ , with cTq = 
1^, Tro-^o-^ = dx6^y, and Tra^ = for j = I, . . .dl - I, provides 
a basis of XCK^). We refer to terms of the form erf' V' (/ > 1) 
as of the type Ai, terms of the form erf' crf^ V^^^ (/, j > 1) 
as of the type A1A2, and so on. In the Appendix, we prove that a 
matrix satisfies condition ^ iff it contains the terms listed in this 
table. Each of the terms can allow signalling in at most one direction 
and can be realized in a situation in which either Bob's actions are 
not in the causal past of Alice's (B ^ A) or vice versa (A ^ B). 
The most general unidirectional process is a quantum channel with 
memory. Measurements of bipartite states that lead to non- signalling 
probabilities can be realized in both situations. The most general 
process matrix can contain terms from both rows and may not be 
decomposable into a mixture of quantum channels from Alice to Bob 
and from Bob to Alice. 

1 - ^ is represented by a process matrix of the form 

We will call processes of this kind causally separable (note 
that the decomposition ([6]) need not be unique since non- 
signalling processes can be included either in W^^^ or in 
W^^^). They represent the most general bipartite quantum 
processes for which the local experiments are performed in 
closed laboratories embedded in a definite causal structure. 
In particular, they generate the most general quantum cor- 
relations between measurements that take place at definite 
(though possibly unknown) instants of time. Clearly, accord- 
ing to the argument presented earlier, causally separable pro- 
cesses cannot be used by Alice an Bob to violate the causal 
inequality ([2]). 

In the Appendix, we provide a complete characterization of 
process matrices via the terms allowed in their expansion in a 
Hilbert-Schmidt basis, which we relate to the possible direc- 
tions of signalling they allow (see Fig. [Sj. We also provide 
possible interpretations of the terms that are not allowed in a 
process matrix (see Fig.|4]and Fig. [5]). 

A causally nonseparable process 

The question whether all local quantum experiments can 
be embedded in a global causal structure corresponds to the 
question whether all process matrices are causally separable. 
Note that this is not a question about entanglement: all possi- 
ble entangled states, and more generally all quantum circuits, 
correspond to matrices of the form W^^^ or W^^^, while the 
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Fig. 4: Terms not appearing in a process matrix. These terms 
are not compatible with local quantum mechanics because they yield 
non-unit probabilities for some completely positive trace-preserving 
maps. A possible interpretation of these terms within our frame- 
work is that they correspond to statistical sub-ensembles of possible 
processes. For example, terms of the type A2 can be understood as 
postselection. One specific case is when a system enters a laboratory 
in a maximally mixed state, is subject to the map M and, after going 
out of the laboratory, is measured to be in some state \if/). The corre- 



sponding probability is given by Tr |^|<A>(<A|A1(^)J, generated in our 

formalism by W^^^^ = -j Notably, correlations of the 

type A1A2 have been exploited in models for describing closed time- 
like curves (43]|45l. The pictures are only suggestive of the possible 
interpretations. 



is remarkable, since if Alice and Bob perform their experi- 
ments inside laboratories that they believe are isolated from 
the outside world for the duration of their operations (e.g. by 
walls made of impenetrable material), and if they believe that 
they are able to freely choose the bits a, b, and b' (e.g. by 
tossing a coin), they will have to conclude that the events in 
their experiment do not take place in a causal sequence. In- 
deed, the framework only assumes that the local operations 
from the input to the output system of each party are correctly 
described by quantum mechanics, and it is compatible with 
any physical situation in which one would have all the rea- 
sons to believe that each party's operations are freely chosen 
in a closed laboratory. 

Interestingly, both the classical bound ^ and the quantum 
violation ([8]) match the corresponding numbers in the CHSH- 
Bell inequality |31 1, which strongly resembles inequality ([2]). 
However, the physical situations to which these inequalities 
correspond is very different: Bell inequalities can be vio- 
lated in space-like separated laboratories, while ([8]) cannot be 
achieved neither with space-like nor with time-like separated 
laboratories. It is an open question whether ^ is the maximal 
possible violation allowed by quantum mechanics. 



non- separable processes we are looking for cannot be written 
as quantum circuits or even as probabilistic mixtures of dif- 
ferent circuits. Surprisingly, an example of such a kind exists. 
Consider the process matrix 



y^AiA2BiB2 



1 



JA1A2B1B2 



1 



(7) 



where Ai, A2, Bi, and B2 are two-level systems (e.g. the spin 
degrees of freedom of a spin-^ particle) and cTx and cr^ are 
the Pauli spin matrices. It can be verified straightforwardly 
that conditions Q and ^ are satisfied, hence ([7]) is a valid 
bipartite process. Having such a resource, Alice and Bob can 
play the game described above and exceed the bound on the 
probability of success ^ imposed by causal order. Indeed, if 
Bob measures in the z basis and detects one of the states |z±), 
the corresponding CJ operator contains the factor |z±)<z±|^^ 
Inserting this, together with Eq. ([7]), into the expression ^ 
for the probabilities, the term containing erf ^ in the process 
matrix is annihilated and what remains corresponds to a noisy 
channel from Alice to Bob. If Alice encodes her bit in the z 
basis with the CJ operator \z±}{z±\^^, this channel allows Bob 
to guess Alice's bit with probability P(y = a) = If, on 

the other hand. Bob measures in the x basis, Eq. ([T]) is reduced 
to a similar noisy channel from Bob to Alice. Bob is thus able 
to activate a channel in the desired direction by choosing the 
measurement basis (see Appendix for a detailed calculation 
and analysis of the protocol). In this way they can achieve 



2+ V2 



^4' 



(8) 



which proves that ^ is not causally separable. We see that, 
depending on his choice. Bob can eflTectively end up 'before' 
or 'after' Alice, each possibility with a probability V2/2. This 



Classical processes are causally separable 

It is not difl&cult to see that if the operations of the local 
parties are classical, they can always be understood as tak- 
ing place in a global causal structure. Classical operations 
can be described by transition matrices M^.'^^'^^^ = P(A2, jl^ii), 
where P(A2, j\Ai) is the conditional probability that the mea- 
surement outcome j is observed and the classical output state 
A2 is prepared given that the input state is Ai. They can be 
expressed in the quantum formalism as CP maps diagonal in 
a fixed ('pointer') basis, and the corresponding CJ operators 
are Mj = Y.a,a, Mf^'\Ai){Ai\^' |^2)<^2l^^ Thus, in order 
to express arbitrary bipartite probabilities of classical maps, it 
is suflftcient to consider process matrices which are diagonal 
in the pointer basis. In the Appendix, we provide a detailed 
proof that all such processes are causally separable. 



Discussion 

We have seen that by relaxing the assumption of definite 
global causal order and requiring that the standard quantum 
formalism holds only locally, we obtain the possibility for 
global causal relations that are not included in the usual for- 
mulation of quantum mechanics. The latter is reminiscent of 
the situation in general relativity, where by requiring that lo- 
cally the geometry is that of flat Minkowski space-time, one 
obtains the possibility of having more general, curved space- 
times. 

The natural question is whether "non-causal" quantum cor- 
relations of the kind described by our formalism can be found 
in nature. One can speculate that they may exist in unprobed 
physical regimes, such as, for example, those in which quan- 
tum mechanics and general relativity become relevant. In- 
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deed, our result that classical theories can always be under- 
stood in terms of a global causal structure suggests the pos- 
sibility that the observed causal order of space-time might 
not be a fundamental property of nature but rather emerge 
from a more fundamental theory |32-34| in a quantum- to- 
classical transition due to, for example, decoherence |[35ll or 
coarse-grained measurements f36l. Once a causal structure 
is present, it is possible to derive relativistic space-time from 
it under appropriate conditions (371 Ell • Furthermore, since 
the conformal space-time metric is a description of the causal 
relation between space-time points |l39j|40l, one can expect 
that an extension of general relativity to the quantum domain 
would involve situations where different causal orders could 
coexist "in superposition". The formalism we presented may 
offer a natural route in this direction: based only on the as- 
sumption that quantum mechanics is valid locally, it yields 
causal relations that cannot be understood as arising from a 
definite, underlying order. 

It is also worth noting that exotic causal structures already 
appear in the classical theory of general relativity. For exam- 
ple, there exist solutions to the Einstein equation containing 
closed time-like curves (CTCs) |41 ]. In this context, it should 
be noted that any process matrix W in our framework can be 
interpreted as a CPTP map from the outputs, A2, B2, of the 
parties, to their inputs, Ai, Bi.ln other words, any process can 
be thought of as having the form of a CTC, where information 
is sent back in time through a noisy channel (see also Fig.[T])). 
The existence of processes that do not describe definite causal 
order is therefore not incompatible with general relativity in 
principle. It is sometimes argued that CTCs should not exist 
since they generate logical paradoxes, such as an agent go- 
ing back in time and killing his grandfather. The possible so- 
lutions that have been proposed |42-47|, in which quantum 
mechanics and CTCs might coexist, involve non-linear exten- 
sions of quantum theory that deviate from quantum mechanics 
already at the level of local experiments. Our framework, on 
the other hand, is by construction linear and in agreement with 
local quantum mechanics, and yet paradoxes are avoided, in 
accordance with the Novikov principle BSl . due to the noise 
in the evolution 'backward in time'. 

Finally we remark that instances of indefinite causal orders 
may also emerge in situations closer to possible laboratory 
implementations. As already noted, our formalism describes 
more general correlations than those that can be realized with 
a quantum circuit, that is, as a sequence of quantum gates. 
Recently, a new model of quantum computation which goes 
beyond the causal paradigm of quantum circuits by using su- 
perpositions of the 'wires' connecting different gates was pro- 
posed |49 |. This possibility may allow breaking assumption 
CS that events are localized in a causal structure. Since the 
instant when a system enters a device depends on how the de- 
vice is wired with the rest of the computer's architecture, su- 
perpositions of wires may allow creating situations in which 
events are not localized in time (similarly to the way in which 
a quantum particle may not be localized in space). While it is 
an open question whether violating the causal inequality ^ 
can be achieved by similar means, the present work suggests 
that new quantum resources for information processing might 



be available — beyond entanglement, quantum memories, and 
even 'superpositions of wires' — and the formalism introduced 
provides a natural framework for exploring them. 
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Appendix 

A. Formal derivation of the causal inequality 

A causal structure (for instance, space-time) is a set of event 
locations equipped with a partial ordering relation < that de- 
fines the possible causal relations between events at these lo- 
cations. If A and B are two such locations, A < B reads "A is 
in the causal past of B'\ or equivalently, "B is in the causal 
future of A" (e.g. if A and B are space-time points, A < B 
corresponds to A being in the past light cone of B). Opera- 
tionally, iiA<B, an agent at A can signal to an agent at B by 
encoding information in events at A that get correlated with 
events at B which the other agent can observe. (Formally, sig- 
nalling from A to 5 is the existence of statistical correlations 
between a random variable at A which can be chosen freely, 
and another random variable at B. By definition, a freely cho- 
sen variable is one that can be correlated only with variables 
in its causal future. Note that a freely chosen variable is an 
idealization since the result of a coin toss or any other candi- 
date for a freely chosen variable may be correlated with initial 
conditions in the past or with space-like separated events, but 
these correlations are ignored as not relevant to the variables 
of interest.) The fact that the relation ^ is a partial order means 
that it satisfies the following conditions: \) A < A (reflexiv- 
ity); 2) '\f A < B din&B < C, then A < C (transitivity); and 
3) if A ^ 5 and 5 ^ A, then A = B (antisymmetry). The last 
condition says that if A and B are two different locations, there 
can either be signalling from A to B, or vice versa, but no sig- 
nalling in both directions is possible (i.e. there are no causal 
loops). If A is not in the causal past of B, we will write A^B. 
Note that in a causal structure both A^B and B ^ A may hold 
(as in the case when A and B are space-like separated), and at 
least one of the two must hold for A B. We will denote the 
situation where both A ^B din&B ^ A hold by A^^B. 

Since every event specifies an event location, we will use 
the same notation directly for events. For instance, if X and Y 
are two events such that the location of X is in the causal past 
of the location of 7, we will write X <Y (similarly for ^ and 
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The main events in our communication task are the systems 
entering AHce's and Bob's laboratories, which we will denote 
by Ai and Bi, respectively, and the parties producing the bits 
a, b, b\ X, and which we will denote by the same letters as 
the corresponding bits. The fact that Alice generates the bit a 
and produces her guess x after the system enters her labora- 
tory means that Ai < a,y. Similarly, we have Bi < b\ b,y. 

The assumptions behind the causal inequality are: 

Causal structure (CS) — The events Ai, Bi, a, b, b\ x,y are 
localized in a causal structure. 

Free choice (FC) — Each of the bits a, b, and b' can only be 
correlated with events in its causal future (this concerns only 
events relevant to the task). We assume also that each of them 
takes values or 1 with probability 1/2. 

Closed laboratories (CL) — x can be correlated with b only 
ifb<Ai, and y can be correlated with a only ifa<Bi. 

We want to show that these assumptions imply 

Psucc = \p(x = b\b' = 0) + ^p(y = aW = 1) ^ ^ (9) 

for the success probability that Alice and Bob can achieve in 
their task. 

First, notice that assumption FC implies that the bits a, b, 
and b' are independent of each other (CS is assumed through- 
out). Indeed, there are two general ways in which the three 
bits could be correlated — two of them are correlated with each 
other while the third one is independent, or each of them is 
correlated with the other two. In the first case, the free-choice 



= ^p(x = bW = 
+ ^p(y = a\b' = 



If Ai ^ ^1 (which implies Bi ^ Ai), from the transitivity of 
partial order it follows that Ai < b (and thus b ^ Ai). From 
assumption CL, x can only be correlated with b if b is in the 
causal past of Ai, thus p(b\x;Ai < Bi) = p(b\Ai < Bi) = \ 
[the last equality follows from the independence of b from the 
causal relations between Ai and Bi, together with assumption 
FC] . Using also that b and b' are independent, we thus obtain 
p(x = bW = 0;Ai < Bi) = p(b = 0',x = 0\b' = 0;Ai < 
Bi) + p(b=l,x= W = 0;Ai < Bi) = p(b = 0\x = 0;b' = 
0;Ai < Bi)p(x = 0\b' = 0;Ai < Bi) -h p(b = l\x = Ub' = 



assumption implies that the two correlated bits would have to 
be in each other's causal pasts, which is impossible. In the 
second case, each of the bits would have to be in the causal 
past of the other two, which is again impossible. Hence, the 
bits are uncorrected. 

Next, consider the following three possibilities that can be 
realized in a causal structure (CS is assumed throughout): 
Ai < Bi, Bi < Ai, Ai Bi. Since these possibilities are 
mutually exclusive and exhaustive, their probabilities satisfy 
p(Ai < Bi) + p(Bi < Ai) + p(Ai Si- Bi) = 1. From assump- 
tion FC it follows that the bits a, b, and b' are independent of 
which of these possibilities is realized. To see this, consider 
for instance b\ Since Bi < b\ we have that b' must be in- 
dependent of whether Ai takes place in the causal past of Bi 
or not, i.e. p(Ai < Bi\b') = p(Ai < Bi). Similarly, b' must 
be independent of whether Ai takes place in the larger region 
which is a complement of the causal future of , which im- 
plies p(Bi i Ai\b') = p(Bi i Ai). But p(Bi i Ai\b') = 
p(Ai < BiW) + p(Ai Si BiW) = p(Ai < Bi) + p(Ai Si 
Bm. while p(B^ i Ai) = p(Ai < B,) + p(A, Si ^i), 
which implies p(Ax Si Bx\b') = p(Ai Si Bi). Finally, 
since p(Ai < BiW) + p(Ai Si BiW) + p(Bi < AiW) = 
p(Ai < Bi) + p(Ai a Bi) + p(Bi < Ai\b') = 1 = p(Ai < 
Bi) + p(Ai Si Bi) + p(Bi < Ai), we have p(Bi < Ai\b') = 
p(Bi < Ai). An analogous argument shows that a and b are 
also independent of the causal relation between Ai and Bi. 

Using the above, the success probability can be written 



0;Ai < Bi)p(x = W = 0;Ai < Bi) = \p(x = 0\b' = 0;Ai < 
B,)^'^p(x=W = 0;A,<B,)='^. 



If Bi < Ai (which implies Ai ^ Bi), by an analogous ar- 
gument we obtain p(y = a\b' = < Ai) = ^. Finally, if 
Ai M Bu we have both p(y = a\b' = 1; Ai M Bi) = | and 
p(x = bW =^\Ai Si Bi) = \. Substituting this in Eq. ( [T0| , 
we obtain 



Psucc =\P(^ = ^\^' = 0) + ^M}^ = = 1) 

0; Ai < Bi)p(Ai < Bi) + ^-p(x = bW = 0;Bi< Ai)p(Bi < AO + ^p(x = bW = 0; Ai M Bi)p(Ai Si Bi) 
1; Ai < Bi)p(Ai < Bi) + ^p(y = aW = UBi< Ai)p(Bi < AO + ^p(y = aW = 1; Ai H Bi)p(Ai Si Bi) 
= (^^p(x = bW = 0; Ai < Bi) + ^-piy = aW = 1; Ai < ^oJ^Ai < Bi) 
+ i^-p(x = bW = 0;Bi< AO ^\p(y = = UBi< Aoj p(Bi < AO 
+ |^/7(x = bW = 0;Ai M Bi) + ^p(y = aW = l;Ai M ^oj^Ai M Bi). (10) 

I 
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Psucc = 1^ + \p(y = aW = l;Ai < Bi)^p(Ai < Bi) 
+ I^M^ = b\b' = 0;B,< Ai) + i j p(B, < A,) 

< ^-p{Ai < B,) + ^-p{B, < AO + ^-p{A, M Bi) = ^. (11) 
This completes the proof. 

B. Definition of process matrices 

In this section we will derive the linear representation ^ 
as well as the conditions ^ and ^ that a process matrix has 
to satisfy. 

Linearity of probabilities. A quantum instrument 1 26 1 is 
defined as a set {MjVj^^ of CP maps such that M = 
is a CPTP map. Our main assumption is that the description 
of the operations in the individual laboratories is in agreement 
with quantum mechanics. In particular, we derive linearity 
from the quantum mechanical representation of probabilistic 
mixtures and of coarse-graining of operations. Consider first 
an instrument [MjY-^^ defined as the randomization of two 
diff'erent instruments {MjYj^^ and {NjY^j^^, where the first is 
performed with probability p and the second with probabil- 
ity (1 - p). The probability to observe the outcome j is, by 
definition, P{Mj) = pP(Mj) + (1 - p)P(Nj). In quantum me- 
chanics randomization is described as a convex linear com- 
bination, Mj = pMj + (1 - p)Nj. We can then conclude 
that the probability must respect linear convex combinations: 
p[pMj + (1 - p)Ni) = pP{Mj) + (1 - p)P{Mj). Consider 
then the coarse-graining of an instrument [MjTj^y This is 
realized when two or more outcomes, for example those cor- 
responding to the labels j = n - \ and j = n, are treated as a 
single one. In the resulting instrument {A1j}"^J all non coarse- 
grained outcomes correspond to the original CP maps Mj = 
Mj for j = 1, ... ^ - 2, while the probability of the coarse- 
grained outcome is given by P(Mn-i) = P(Mn-i)-\-P(Mn). In 
quantum mechanics, the CP map corresponding to the coarse 
graining of two outcomes is represented by the sum of the 
respective CP maps, Mn-i = Mn-i + At„, from which it fol- 
lows that P (Mn-i + Mn) = P (Mn-i) + P {Mn)- Randomiza- 
tion and coarse graining together impose linearity. The argu- 
ment can be repeated for two (or more) parties, yielding the 
conclusion that all bipartite probabilities compatible with a 
local quantum mechanical description are bilinear functions, 
P(Mf,M^) = co[Mf,Mf) e [0,1], of the local CP and 
trace-nonincreasing maps Mf, Alf . 

Thanks to the CJ isomorphism, it is possible to represent bi- 
linear functions of CP maps as bilinear functions of matrices: 
oj ^ (I) : £(9^^' X £(9(^' ^K^^) ^ R. In general, 
multilinear functions on a set of vector spaces V ^ x x . . . are 



isomorphic to linear functions on V^^V^®. . . , hence the prob- 
abilities can be written as linear functions on X(*H^^ *H^^ 
TY^i (g) T^^O- Using the Hilbert-Schmidt scalar product, we 
can identify each real linear function with an element of the 
same space, oj ^ w^'^^^'^^ e £(9{^^ ^K^^ ^ <j^b, ^ ^B2^^^ 
arriving at the representation ([3j. 

Nonnegativity and normalization of probabilities. The re- 
quirement that the probabilities are non-negative for any pair 
of CP maps and imposes the restriction that W is pos- 
itive on pure tensors (POPT) 1 50] with respect to the partition 
A1A2 - B1B2. These are matrices such that 

^ ^ (12) 

The condition has to be imposed for arbitrary positive 
semidefinite matrices M^^^^ and M^^^^ because these are the 
CJ matrices of CP maps. 

We additionally assume that the parties can share arbitrary 
(possibly entangled) ancillary states independent of the pro- 
cess, and use them in their local operations. The latter means 
that each party can extend the input space of her/his opera- 
tions to the ancillas, which we denote by A[ and B[ for Al- 
ice and Bob, respectively, and apply arbitrary quantum op- 
erations with CP maps : XCK^i ^ £('H^'), 
: £{'H^'^ ^ £{'H^^). (One can similarly ex- 
tend the output systems, but this is not necessary for our ar- 
gument.) The assumption that the ancillary systems contain a 
joint quantum state independent of the process means that if 
separate operations are applied on the ancillas and the original 
systems, the joint probability distribution for the outcomes is 
a product of two distributions — one for the outcomes on the 
ancillas, which is the same as one arising from a measurement 
on a quantum state p^'i^'i , and another one for the outcomes on 
the original systems, which is given by Eq. ([3]) with the orig- 
inal VK^i^2fiifi2 These requirements imply that the extended 
process matrix is given by w^'^^'^^^'^^'^^ = p^i^i iy^i^2fiifi2 
If we then require that the probabilities for extended opera- 
tions are non-negative, one has 

Tr [p^i^i W^i^2fiifi2 (^M^[MA2 ^ j^B[B,B2)^-^ ^ ^^3) 

^^J^A[A,A2^J^B[BlB2^^A[B[ 

It was shown 1 50 1 that condition ([TSj is satisfied if and only if 
^^AiA2BiB2 positive semidefinite (a class strictly smaller than 
POPT), which is condition (|4]). 

Additionally, probabilities must be normalized: 1 = 
Y^ij OJ (Mf , M^) = OJ (2, Mf , M^), which means 

oj[Af^,M^) = 1, VCPTPM^, M^. (14) 

Condition ^ can be deduced from Eq. ([14]) simply by notic- 
ing that for a CPTP map M the corresponding CJ matrix sat- 
isfies the condition Tr^^M^^^^ = Tr^^ (J Al(|0+)(0+|))'^ = 
[TrA2 (1 M(\(f)^}{(f)^\))f = 1^^ . To see that this is also a suf- 
ficient condition for a map to be trace-preserving, it is enough 
to consider the inverse direction of the CJ isomorphism, 

M(p^i) := (tta, [p^^M^^^^])^ . (15) 



9 



C. Characterization of process matrices 

Here we derive necessary and sufficient conditions for a 
matrix w^^^^^^^^ to satisfy Eq. ^ and Eq. ([Sj in terms of an 
expansion of the matrix in a Hilbert- Schmidt basis. A Hilbert- 

Schmidt basis of Jl^H^) is given by a set of matrices {cr^f^^^Q , 
with cr^ = tx, T^^cr^CTy = dxS^y, and Trcrj = for j = 
1, ... 4 - 1. A general element of £CH^' ^K^^ 
can be expressed as 



(16) 



livAy 



(we omit tensor products and identity matrices whenever there 
is no risk of confusion). Since a process matrix has to be 
Hermitian, we consider only the cases 



(17) 



We will refer to terms of the form cr^^ V^^^ (^ > 1) as 
of the type Ai, terms such as cr^^ cr^^ f^^^ (/, 7 > 1) as 
of the type A1A2, and so on. The properties of a process ma- 
trix can be analysed with respect to the terms it contains. For 
example, terms of the type AiBi produce non-signalling cor- 
relations between the measurements, terms such as A25i cor- 
relate Alice's outputs with Bob's inputs, yielding signalling 
from Alice to Bob, etc., as illustrated in FigJS] Note that not 
all terms are compatible with the condition ([5f7 We will prove 
that a matrix W satisfies condition ([5]) if and only if it only 
contains the terms listed in Fig. [3] 

The CJ matrix of a local operation can be similarly written 
M^i^2 ^ Y,^^ ^^^^x,^X2^ ^ condition Tr^^M^^^^ = 



is equivalent to the requirement roo = Ho = for 
/ > 0. Thus CJ matrices corresponding to CPTP maps have 
the form 



1 



ij>0 



i>0 

aiJij e R. 



Let us consider first the case of a single party, say, Alice. 
Since the set of matrices M^^^^ > is a substantial set, condi- 
tion ([5]) can be equivalently imposed on arbitrary matrices of 
the form ([18]) and, for a single party, it can be rewritten as 



1 T 

"A2 



/>0 ij>0 

V ai, tij e R. 



1, 



Using an expansion of the process matrix in the same basis in 
a similar way, iy^i^2 - ^^^y w^ycr^ ^ w^y e R, the above 
condition becomes 



Woo + ^ y^OiCli + ^ Wijtij 
i>0 ij>0 

V ai, tij G R, 



= 1, 



and one obtains woo = ^ , wo/ = Wij = for /, j > 0. Thus the 
most general process matrix observed by a single party has the 



form 



Vi e R, > 



(19) 



which can be recognized as a state. This result — that all prob- 
abilities a single agent can observe are described by quantum 
states — is an extension of Gleason's theorem from POVMs 
ISTl l52l to CP maps (note that here the linear structure of 
quantum operations is assumed, while in Gleason's theorem 
for POVMs it is derived from diff'erent hypotheses. However, 
by a similar argument one could derive linearity for CP maps 
too). 

Let us now consider a bipartite process matrix, \y^i^2fiifi2 - 
l^/uvAy "^fivAyCrfi^cry^o-fcry^ w^y;iy G R. We have to impose ([5j 
for arbitrary matrices M^^^^, M^^^^ of the form ([18]). First, if 
we fix M^i^2 = we obtain 

dBj ' 



= 1 



>^0000 + ^ >^0/00^/ + ^ >^/jOO^/j 
< />0 ij>0 

V ai, tij e R, 

which imposes woooo = and wo/00 = w,/oo = for /, j > 

0. Similarly, by fixing M"^^^^ = we can derive wooo/ = 

^ooij = for /, j > 0. Finally, imposing ([5]) for arbitrary 



A1A2 



1 



i>0 



ij>0 



B1B2 



^2 V k>o ki>o ) 



(18) we obtain 



^ yvoiokaibk + ^ woikiaiSki 

ik>0 ikl>0 

+ ^ ^ijOktijbk + ^ WijkltijSkl = 0, 
ijk>0 ijkl>0 

V ai,tij,bk. Ski ^ 

from which we conclude that the most general matrix that sat- 
isfies (jSj has the form 



y^AiA2BiB2 
^B<A 

<B 



dA,dB, ^ ' 

i7>0 iik>^ 



ij>0 ijk>0 
/>0 />0 ij>0 

where Qj, dijk, eij, fyk, gij, v,-, Xi e R. 
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This form, together with the condition ^^1^2^152 > com- 
pletely characterizes the most general bipartite process ma- 
trix. 



D. Terms not appearing in process matrices 




The not-allowed terms are listed in Fig.|4] along with pos- 
sible interpretations. Particularly interesting are the cases in- 
volving terms of the type A1A2. These would correlate Al- 
ice's output with her input and not give unit probabilities 
for some CPTP maps that she can choose to perform. This 
kind of correlations resemble a 'backward in time' transmis- 
sion of information: one can imagine that they can be gen- 
erated by a quantum channel 'in the inverse order', from 
the output A2 to the input Ai. It is worth noting that a re- 
cently proposed model of closed time-like curves ||43l |45l 
can be expressed precisely in this way. Using our terminol- 
ogy, such a model considers an agent receiving two quantum 
systems in her laboratory: a chronology-respecting system A 
and a second system A' which, after leaving the laboratory, 
is sent back in time to the laboratory's entrance (see Fig. [5]). 
This can be described by the process matrix W^^^'^^^^'^ = 
0-^1 t\(f)+}{(f)+fi^2 (g) u), where o^' is the state 
of the chronology -respecting system when it enters the labo- 
ratory and ]l|0^)(0^|^'i'^2 ^7^ (g) 1^ is a process matrix cor- 
responding to a unitary U from A'^ to A\ , describing the evo- 
lution back in time of the chronology-violating system. (The 
labels Ai, A'^ represent the two systems entering the labora- 
tory, while A2, A2 represent the systems going out. Note that 
here the two systems belong to the same laboratory and they 
can undergo any joint operation.) In this model, probabilities 
have to be renormalized in order to sum up to one, which in- 
troduces a non-linearity that violates our original assumptions 
(in particular, as opposed to quantum mechanics, probabilities 
are contextual in this model, since it is necessary to specify the 
events that did not occur in order to perform the renormaliza- 
tion step). The same can be said for Deutsch's model of closed 
time-like curves |42|, which is also non-linear (although it 
uses a different mechanism to obtain well-defined probabili- 
ties) and thus violates our premise that ordinary quantum me- 
chanics holds locally in each laboratory. 



E. Violation of the causal inequality 

The process described by Eq. ^ can be exploited for the 
task described above in the following way. Alice always mea- 
sures the incoming qubit in the z basis, assigning the value 
X = to the outcome \z+) and x = 1 to She then repre- 
pares the qubit, encoding a in the same basis, and sends it 
away. It is easy to see that the CP map corresponding to the 
detection of a state and repreparation of another state |0) 
has CJ matrix \(p){(p\^^. Accordingly, the possible 

operations performed by Alice can be represented compactly 



Fig. 5: Nonlinear model of closed time-like curve. In the model of 
closed time-like curves considered in Refs. 14311451 , a chronology- 
respecting system A, initially in a state cr, interacts with a second 
system, A' , which travels back in time according to a unitary U. This 
model can be represented in our formalism by an 'unphysical' pro- 
cess matrix, i.e. one for which probabilities do not sum up to one. 



by the CJ matrix 



e'^'Kx^a) = - [H + (-l)V,f ^ [1 + (-l)V,] 



A2 



(20) 



Bob adopts the following protocol. If he wants to read Alice's 
bit (h' = 1), he measures the incoming qubit in the z basis and 
assigns y = 0, y = 1 to the outcomes \z-}, respectively 
(the repreparation is unimportant in this case). If he wants to 
send his bit (b' = 0), he measures in the x basis and, if the 
outcome is he encodes b in the z basis of the outgoing 
qubit as ^ 1 \z-}, while, if the outcome is he 
encodes it as ^ 1 ^ The CJ matrix representing 
Bob's CP map is 



') = b'T]f'^'(y, b) + {V e l)r]^'^\y, b\ 



(21) 
(22) 

4^^'(y,b) = i [1 + (-l)V,]^^ [1 + (-l)^^V,]^' , 

(23) 

where p^^ arbitrary state prepared when Z?' = 1 (with 

Trp^2 - \^ 2ind © denotes the sum modulo 2. Note that in Eq. 
( [23] ) Bob's assignment \x+) ^ y = 0, \x-) ^ y = 1 for the 
outcome of his measurement is arbitrary since for b^ = he is 
not trying to correlate y with a. 

The probabilities for diff'erent possible outcomes, 
when the described protocol is applied to the pro- 
cess (|7]), are given, according to ([3]), by P{xy\abb') = 
Tr ^w^iMBiB2 (^^MA2^-^^ a)r]^'^^(y, b, b'))\. In order to 
calculate the success probability, we need as interme- 
diate steps P(y\ab,b' = 1) = Y^x Pi^yl^^^^' = 1) and 
P(x\ab,b' = 0) = Y.yP(xy\ab,b' = 0). Notice that when 
the outcome of one party is ignored, it is always possible to 
identify a specific state in which the other party receives the 
qubit. For example, to average out Alice's outcomes one has 
to calculate 



^ Tr [w^^^^^^^^ {^^^^^(x, a)rj^^^^(y, b, b'))] 



= TvB,B2\ri^'^'(y.b,b')TTA,A2 



W 



A1A2B1B2 
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The process observed by Bob is therefore described by the 
reduced matrix 



(24) 



The matrix 2x ^^^^^(x, ci) represents the CPTP map performed 
by AHce when the outcomes of her measurement are ig- 
nored (the expHcit dependence on a accounts for the pos- 
sibility of signalHng). Using ( [2Q| we find Y^x^^^^^^^^^) - 



[1 + (-If cr^]^% which, plugged into Eq. ([24]) together 



with Eq. ([7]), gives 



w^i^V . 1 



]l + (-ir— cr, 
V2 



\B2 



(25) 



When this is measured with the map ([22]), we find 



P(y\ab,b' - 1) = Tr 



j]f'^\y,b)W^'^\a) 



1 

2 

2+ V2 



1 -h 



V2 



(26) 



from which we obtain P(y - a\h' = 1) = 

Consider now the case when h' - ^. When Bob's out- 
comes are ignored, he performs the CPTP map described by 
l]yr]^'^'(y,b) =^11 + (-l)Vfvf^]. From this we can cal- 
culate, as in the previous case, the effective state received by 
Alice, which is 



W^'^\b,b' = 0). 



1 



1 

^2 



1^2 



(27) 



from which we find P(x = b\b' = 0) = In conclu- 

sion, the protocol described yields the probability of success 
([8]), which proves that the process in Eq. ([?]) is not causally 
separable. 



F. Casual order in the classical limit 

Let us now show that in the classical limit all correlations 
are causally ordered. Classical operations can be described 
by transition matrices M^^^^ = P(k, where P(k, is the 
conditional probability that the measurement outcome j is ob- 
served and the classical output state k is prepared given that 
the input state is /. They can be expressed in the quantum 
formalism as CP maps diagonal in a fixed ('pointer') basis, 
and the corresponding CJ matrices are Mj = Yjki 
I^X^I'^^ In order to express arbitrary bipartite probabilities of 
classical operations, it is sufl&cient to consider process matri- 
ces of the standard form 



(28) 



where cr^^^ and cr^^^ are diagonal in the pointer basis. Prob- 
abilities are still given by 

P (Mf , M^) = Tr [iy^i^2fiifi2 (^j^MA2 ^ M^'^')] . (29) 



We will show that any such diagonal process matrix can be 
written in the form 

^^A,A2BlB2 ^ I /pAiA2fii _^^MB,B2\ ^ ^3q>^ 

dAidB, ^ ^' 

where p^i^2fii and p^i^i^2 ^^q. positive semidefinite matrices. 
This is sufficient to conclude that w^^^^^^^^ is causally separa- 
ble. Indeed, if ^^"^1^25152 could be written in the form ( [3Q| , we 
know that p^i^2fii would not contain Hilbert- Schmidt terms 
of the types A1A2 or A2 (which are not allowed in a pro- 
cess matrix), since by assumption these terms are not part of 
W'^''^^^'^^ Therefore, the matrix 



r,AiA2Bi 



(31) 



which is positive semidefinite, has trace (iA2^fi2' contains 
only terms of the allowed types, would be a valid process ma- 
trix with no signalling from Bio A. Similarly, 



nAiBiB2 



Trp^ 



B1B2 



(32) 



would be a valid process matrix with no signalling from A 
to B. The whole process matrix could then be written in the 
causally separable form 



^MA2B,B2 ^ ^^BM + (1 _ q)W^iB^ 



where 



q = 



A1A2B1 



dAydA2dBidB2 



(33) 



(34) 



Note that < ^ < 1 since p^i^2fii ^nd p^i^i^2 Eq. ([30]) are 
positive semidefinite and TvW^^^^^^^^ = JA2^fi2- 

To prove Eq. ( [3Q| , we will construct p^i^2fii ^^d p^i^i^2 
from the general form in Eq. ( [28] ). Let the minimum eigen- 
value of cr^^^ -h cr^^^ be m. Since iy^i^2fiifi2 positive 
semidefinite and cr^^^ + cr^^^ is traceless, we have m e 
[-1,0]. Define the matrices 



f^iBiB2 ^ ^iB^ 

The full process matrix can then be written 



(35) 
(36) 



W 



MA2B,B2 ^ ((1 + + /^^^^^ + K^^'^'A , (37) 

dA^dB ^ ' 



where k^^^^^^ + ^^^^1^2 positive semidefinite. 

We are now going to modify A^i^2fii ^j^^j j^iBiB2 adding 
matrices of the form k^^^^ to k^^^^^^ and subtracting them 
from /c^i^i^2 (therefore leaving /c^i^2fii + ^iBiB2 unchanged), 
until we transform both a^i^25i ^^d a^i^i^2 Eq. ([37]) into 
positive semidefinite matrices. 

Denote the pointer basis of system X by |/)^, / = 1, J^, 
X = Ai,A2,Bi,B2. All matrices we consider are diago- 
nal in the basis {\i)'^'\j)'^Hk)^'\l)^^}. Let mi(ij,kj) denote 
the eigenvalues of A^i^2fii corresponding to the eigenvectors 
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\i)^'\j)^^\k)^'\l)^\ and let m2(ij,kj) be the eigenvalues of 
^1^1 ^2 corresponding to the same vectors. For every / and 
k, we do the following. Define 

mi(/, k) = minmiO', k, /), (38) 

m2(i, k) = min m2(/, j, k, I). (39) 

Note that mi(i,j,kj) do not depend on / since a^i^2Si ^^ts 
trivially on B2, and similarly m2(i,j,kj) do not depend on 
j. This means that for given / and k, the minimum of the 
eigenvalues of k^^^^^^ + /(^i^i^2 eigenvectors of the type 

\i}^'\j}^^\k}^'\l}^^ is equal to mi(i,k) + m^Uk). But by con- 
struction /(^i^2fii + /^i^i^2 |g positive semidefinite, so we have 

fhiiUk) + miUk) > 0. (40) 

Now, if both mi(i,k) and m2(i,k)} are non-negative, we will 
not modify k^^^^^^ and /(^i^i^2 However, if one of these num- 
bers is negative, say mi(/, k) <0 (both cannot be negative due 
to (|40|), we will add the term -mi(/, k)\i}{i\'^' 1^^ I^X^I^^ 
11^2 to A^i^2fii and subtract the same term from /^i^i^2 Af- 
ter this step, the modified A^i^2fii g^^;.]^ ^j^^^ ^j^g eigenval- 
ues mi(/, J, ^, /) have been changed to mi(/, 7, ^, /) - mi(/, /:) > 



mi(/, k)-mi(i, k) = 0, i.e. a^i^2Si ^qq^ ^ot have any more neg- 
ative eigenvalues mi(/, j, k, I) for the given / and k. The same 
holds for /(^i^i^2 since the eigenvalues m2(/, 7, ^, /) change to 
^liL j, k, /) -h mi(/, /:) > m2(i, k) -h mi(/, ^) > 0. In other words, 
the eigenvalues of the modified k^^^^^^ and k^^^^^^ satisfy 

mi (/, 7, k, /), m2(/, 7, ^, /) > 0, V 7, /. (41) 

By performing this procedure for all / and k, we eventu- 
ally transform A^i^2fii ^iBiB2 ^^^^ matrices all of whose 
eigenvalues are non-negative. Denote the resultant positive 
semidefinite matrices by ^i^2fii ^iBiB2 
the term (1 -hm)! in Eq. ( [37] ) for instance to ^i^25i (recall that 
m G [-1,0]), defining the positive semidefinite matrices 

pA,A2B, = Q + + ^lA^fii ^ (^42) 

pMB,B,^^^,B,B,^ (43) 

We thus arrive at the desired form ( [3Q| which implies ( [33] ) as 
argued above. 
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